On the Asymptotic Stability of De-Sitter Spacetime: a 
non-linear perturbative approach 

o 

O Filipe C. Mena* 

<N 

^ Dep. Matemdtica, Universidade do Minho, 4710-057 Braga, Portugal 

I— i 

^ June 15, 2010 

Oh 
I 

^ Abstract 

We derive evolution and constraint equations for second order perturbations of 
^ flat dust homogeneous and isotropic solutions to the Einstein field equations using all 

scalar, vector and tensor perturbation modes. We show that the perturbations decay 
I asymptotically in time and that the solutions converge to the De-Sitter solution. By 

J>. induction, this result is valid for perturbations of arbitrary order. This is in agreement 

CD with the cosmic no-hair conjecture of Gibbons and Hawking. 
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1 Introduction 



X 

a long standing problem in Mathematical General Relativity is the proof of the Gibbons- 

^ Hawking cosmological attractor conjecture, best known as the cosmic no-hair conjecture. 

This conjecture is the statement that generic expanding cosmological solutions to the Ein- 
stein Field Equations (EFEs) with a positive cosmological constant A tend asymptotically 
in time to the De-Sitter solution [8]. Wald provided a proof of this conjecture within the 
class of spatially homogeneous solutions [16] while other authors, including Stein-Schabes 
and collaborators [1, 9], have used particular classes of inhomogeneous exact solutions 
with symmetries. So, it is important to consider more general non-symmetric cases and 
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this can be achieved by using perturbation theory. An interesting result along this line 
was derived in [12] where the author shows, subject to conditions on a symmetric hy- 
perbolic system, the exponential decay of linear and non-linear perturbations around flat 
homogeneous and isotropic cosmologies. Other results, which do not necessarily imply an 
exponential decay, were obtained for linear metric perturbations around homogeneous and 
isotropic backgrounds by [13, 2, 11]. 

A second order perturbative approach was followed by [3] using scalar perturbations 
modes at first order which then coupled to source second order scalar and tensor pertur- 
bations. In particular, it has been shown that those first and second order perturbations 
decay and that the respective perturbed solution to the EFEs approaches the De-Sitter 
solution asymptotically in time [3]. 

In this note, we generalize the results of [3] in three ways: firstly, by including vector 
and tensor perturbations at first order (besides the scalar modes of [3]). Secondly, by 
considering vector perturbations at second order, so that all scalar, vector and tensor 
modes are included at both first and second orders. Thirdly, by considering all those 
perturbations at any order. To do that, we shall first derive the second order perturbation 
equations around flat dust homogeneous and isotropic backgrounds by generalising the 
work of [10]. We shall then estimate the asymptotic dynamics of the solutions and show 
that they decay in time. This will enable to prove that the perturbed spacetime asymptotes 
to De-Sitter. Finally, we shall show that, by induction, this result is valid for perturbations 
of any order. 

We use the following indice notation: small greek a, /?,... = 0,1,2,3 and small latin 
o, 6, ... = 1, 2, 3. 

We shall consider solutions to the EFEs such that Ta/s = pUaU^, where p is the 
mass density and Uq, the 4-vclocity of the fluid. We denote by E the space orthogonal 
to u on spacetime {M,g). Following [10], we use comoving euclidean coordinates, say 
{i, x^, x^, x^}, to write the components of tt as = (1,0,0,0) and, in this basis, the 
metric on M can be written as: 

g = -dt^ + habdx"'dx\ (1) 

where hat is a function of all coordinates. The components of the extrinsic curvature on S 
coincide with the expansion tensor (see [6] for definitions about the 1-1-3 formalism) and 
can be written as 

where the dot denotes differentiation with respect to t. In this case, the EFEs give (see 
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also [10]): 



(2) 

(3) 
(4) 
(5) 



where we have denoted the components of the Ricci tensor on S as i?*^ and the Ricci 
scalar on S as i?*. 



2 Perturbations around FLRW flat backgrounds 

In this section, we consider small perturbations of arbitrary order n G IN on a flat 
Priedmann-Lemaitre- Robertson- Walker (FLRW) background spacetime (Mb, ^b) with A > 
and Tap = pUaUp. Following [10], we take the conformally rescaled FLRW flat metric 

= ffi*;] = a\T){-dT^ + Sabdx'^dx''), (6) 

where r is the conformal time defined by dt = adr and a is the scale factor. We consider 
perturbations in the form [10] 
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where V^"^ can be decomposed as 



n) 



(7) 

(8) 
(9) 
(10) 

(11) 



with VWi"^ = 0. Now, for this case, as was shown by Stewart [14] (following works of 
Choquet-Bruhat, Fischer & Marsden [4] and D'Eath [5]): 
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Theorem 1 (Stewart) Any symmetric tensor x^/s defined on three-surfaces S of can 
be decomposed into three parts 

xi? = (VaV, - ^5„,A)x(") + 2V(„4") + (12) 

where A = V"Va, the scalar x defined up to a constant, the vector field is 
divergence-free and tt^^ is both transverse and tracefree (i.e. V"7r^^^ = tt^q^" = 0/ 

There are, however, too many degrees of freedom in the above decomposition, and one 
can set some of the free perturbations functions to zero thus choosing a certain gauge. For 
convenience, we consider the synchronous gauge which, in this case, corresponds to taking 

v„$(") = V^''^ = W^i"^ = 0. 

However, we shall prove our main results in a gauge independent manner. 

We note that, in this framework, the initial perturbations (8)-(10) are assumed to be 
small, in the sense that they result from a Taylor series expansion for each order 0{e"'), 
with e ^ 1 (see e.g. [14, 10]), and are such that the metric can be written in "almost- 
RW" form (holding the inequalities (2.58) of [7]). We also note that the dependence of 
the asymptotic dynamics on more general initial conditions cannot be rigorously studied 
with the present methods. In order to do that, one could use methods from the theory of 
hyperbolic partial differential equations, such as the ones used in [11, 12]. 



3 Evolution of first order perturbations 

The evolution of first order perturbations in this context is a much studied subject (see 
e.g. [3] and references therein). Here we shall just recall the main results which will be 
needed for the next section. 

From the EFEs with the perturbed metric (8)-(10) on FLRW flat backgrounds in the 
synchronous gauge one finds that the evolution equations for the first order perturbations 
are: 



a^ 2 

4^^' + -4^^ = (14) 



+ - i,a'p,c^('^ = (13) 



J 



a 



where the prime denotes differentiation with respect to conformal time. Note that, in the 
synchronous gauge, the only non-zero perturbation variables are (/)(^) , 7r^^\ x^^^ but 
(^^^■^ and ')^^ are related through the momentum constraint equation (4). 
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Figure 1: Plot obtained from numerical integration of equation (15) after Fourier decom- 
position for A = 0.001, po = 0.01 and the Fourier mode q = 1 (see equations (18) and 
(19)). 

It is known that the solutions of (13) increase in time and approach a positive constant 
asymptotically. In turn, the solutions of (14) decrease in time and asymptote to zero, 
while the solutions of (15) oscilate and decay asymptotically to a constant value in time 
(an example of this behaviour is depicted in Figure 1). One way to show this is to use the 
asymptotic form for the scale factor which results from the Friedmann equations (see e.g. 
[3]): 

a{r) = -^l- (16) 

Note that as a — >■ +oo one has r ^ 0. So, we are interested in the asymptotic behaviour 
of the solutions to the perturbation equations as r — > 0. Using (13) and (16) one gets 

,^(i)(t,x) = Ci(x)rr (I ^Ar'/'^ + C2(x)rJ ^Ar'/'^ , (17) 



where = j,, Ci and C2 are C^-functions which depend on the initial data and 

J and Y are the Bessel functions of the first and second kind, respectively. One can also 
find a formal solution to (15) which can be written in terms of an integral over the wave 
vector space 

^li^(^,^)=/ 7r«(g,r)e-'i-dq, (18) 

JO 



with 



1) = ^ab{Q){sin {qr) - cos (gr)) + F„6(g)(cos (gr) + sin (gr)), (19) 
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where Kab and Fab are functions (of the Fourier mode q) which depend on the initial data, 
but are otherwise arbitrary. Finally, from (14) one finds 

= (20) 

where C(x) depends on the initial conditions. So, from (18), (19) and (20) one gets that 
the linear perturbations cp^^'^ , n^fj and Z^^^ tend asymptotically in time to constant values. 
This corresponds to a local asymptotic approach to the De-Sitter solution. In order to see 
this explicitly one can use the Boucher-Gibbons [2] local coordinate change to transform 
the resulting asymptotic form of (7) into the De-Sitter metric (see e.g. [3]). This enables 
to show: 

Proposition 1 For n = 1, the metric (7), satisfying the Einstein Field Equations, ap- 
proaches the De-Sitter solution locally asymptotically in time. 



4 Evolution of second order perturbations and non-linear 
asymptotic stability of De-Sitter 

One can carry over to second order the first order analysis. In this section, we present the 
evolution and constraint equations for second order perturbations which include scalar, 
tensor and vector modes at first and second orders in presence of a non-zero cosmological 
constant {A ^ 0). These equations are derived in the synchronous gauge from (2)- (5) and 
generalise the equations obtained in [10] for the case of A = and zero vector perturbation 

(2) 

modes. For convenience, we do not split xib vector and tensor parts as before 

(recall equation (12)) and we write the perturbations in the 3-metric as g^^ = a^{T)j^ 
and g^^ = a^{T)j^^. 

Raychaudhuri equation 



1 2 
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Energy constraint 



+ 



24 



+ 



7 



3 a 



(^v^7^;^+7r::,-27r;.j + 



, ^(l)<ia ( {l)b _ {l)b\ 3 _ I {l)ab (l)d _ 1 (l)a,d (l)^' 
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(22) 
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(2) 



Momentum constraint 

^'1 



2^(2)' + 1 (V _ ^(l)a<if ^(1) ' "\ (1)' _ 1 (l)«<i (1)' _ 1 (l)a (l)d' ^ ^ 

^^,6 ^ 2^^^ ~ ' I 'Ha i ^7 ,a%d ^l ylad ^'a.dlb 

Evolution equation 

-(^(2)%2^^(2)')n+^(xr+24?^"') 

1 (2)da 1 (2)d,a _ ly2 (2)a _ (l)ad' (1)' _ 1 (l)<i' (l)a' 
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(24) 



+ 27^^^*' (7f i + 7^:.^ ,t - ^Id - 7i^ ;) + 2,(^)J - - 

+ 27«rn 7iP'" - 27«rn 7^^? + 7^'^:;7« '» + 7^!e" (7^^^r + 7^'^^'" - 7^ ^"'^ 



(l)dn ra , ^ {l)dn (l)m ™ , 1 (l)d,m (l)„ ™ 

7 ,m7(i„ 0^+7 ,m7d ,„0 6 + .Id In ,m" 



(2)a 



4b 



where S = {p — p^)/ Pb, the subindex denotes evaluation at an initial time tq and the 
commas denote V derivatives. 

Notice that, while at first order the perturbation equations for the difi'erent modes 
(scalar, vector and tensor) decouple from each other (recall equations (13)-(15)), at second 
order the perturbation equations (and, therefore, the perturbation modes) are coupled. 
So, for example, products of first order scalar-tensor, scalar-scalar or scalar-vector modes 
will, in general, source all second order modes. 
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Notice also that the terms N^^^^ on the r.h.s. of the four equations are of the form 
0(1)0(1), i.e. these terms result from couplings of first order perturbations. So, from 
(17)-(20) we get that all terms N^'^'^ approach constant values asymptotically in time. 
This will be crucial ahead to find, from (21)-(24), that all tt^^-* and ^^^^ also tend 

asymptotically to constants in time (an example of this behaviour is depicted in Figure 
2). _ _ 

In detail, using the asymptotic expressions (16)-(20) in equation (21) we can integrate 
to get the asymptotic estimate 

</.(^)(r,x) = ^+0(r^), (25) 
where L{x) is a constant function in time whose expression can be obtained from the 

(2) 

asymptotic form of . In turn, from integration of (22) and (23) we obtain asymptotic 
conditions in the form 

V„xr"(^'^) = Q/3(x) + 0(r2) 
VaV^X^'^"''(r,x) = G(x) + 0(t2) (26) 

where, again, Q/j and G are constant functions in time. Also, substituting the expressions 
(16)-(20) in equation (24), the latter can formally be integrated to give 

/"OO 

xS(^>x) = x^^lir,k)e-^^^dk, (27) 

with 

X^lir,k) = ^ + Oir^) (28) 

where E is a function of the Fourier modes k which depends on the initial conditions. 
From (25)- (28) it is then clear that both ^^^^ and xfb ^^^'^ to constant values in time 

(2) 

asymptotically. Finally, to see that Za individually tends to constant values we substitute 
x'^^^ = 2V^aZ^f + TT^J^ into equation (24). Given Vtt^^ = 0, then there will be no tt^^ 
terms in (24). Now, since the first order terms (on the r.h.s.) and the terms in (p^"^^ (on 

(2) 

the l.h.s.) are already known to asymptote to constants, then the remaining term in Za 
will also tend to a constant value in time. We can use this in equation (24) to show that 

(2) 

TT^j also tends to constants. 

Therefore, we conclude that all second perturbations 0*^^^ , vr^^'* and Z^^ tend asymp- 
totically in time to constant values. Wc note that the density perturbations can be written 
in terms of the metric perturbations and the expressions that relate these quantities have 
been derived in [10]. In the present case, since the metric perturbations asymptote to 
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Figure 2: Two plots obtained from numerical integration of equation (21) for the evolution 
of second order scalar perturbations (p^"^^ for a = FLRW dust background with A = 
0.001 and po = 0.01. 

constants in time, so do the density perturbations. As before, one can now apply the 
Boucher-Gibbons [2] local coordinate change to explicitly transform the resulting asymp- 
totic form of (7) into the De-Sitter metric. 

One can also see this in a coordinate and gauge invariant way: Using (25)- (28) we can 
show that at second order, asymptotically, 

9^ A, Haf3 ^0, 4A, (29) 

0, E^fs -^0, R* ^ 0, (30) 

where a is the shear scalar, H^p and E^p are the magnetic and electric parts of the 
Weyl tensor, respectively. The asymptotic behaviour (29)-(30), together with the fact 
that = and — >■ 0, where oja^ is the vorticity, caracterizes gauge invariantly the 
approach to the De-Sitter solution (see e.g. [7]). 

This result generalises the result of Proposition 1 and previously known second order 
results [3] by including, at second order, sources of first order vector and tensor perturba- 
tion modes as well as pure second order vector modes: 

Proposition 2 For n = 2, the metric ( 7), satisfying the Einstein Field Equations, ap- 
proaches the De-Sitter solution locally asymptotically in time. 

Consider now nth order perturbations in (7). It is clear that, as before, the sources 
of the EFEs for the nth order perturbations are products of lower order terms such as 
0(n - l)O(l), 0(n - 2)0(2) etc. For example: 

</>('^)" + - JpeaV") = ivf \ (31) 
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where A^^ is a sum of products of lower order terms. By proving that the lower order 
terms tend to constants asymptotically in time, one can proceed as before, order by order, 
and show, by induction, the following result: 

Theorem 2 At each order n, the truncated metric (7), satisfying Einstein Field Equa- 
tions, is convergent and approaches the De-Sitter solution locally asymptotically in time. 
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